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Abstract. A time- reversal invariant topological insulator can be generally defined 
by the effective topological field theory with a quantized 9 coefficient, which can only 
take values of or n. This theory is generally valid for an arbitrarily interacting 
system and the quantization of the 8 invariant can be directly measured experimentally. 
Reduced to the case of a non-interacting system, the 6 invariant can be expressed as an 
integral over the entire three dimensional Brillouin zone. Alternatively, non-interacting 
insulators can be classified by topological invariants defined over discrete time-reversal 
invariant momenta. In this paper, we show the complete equivalence between the 
integral and the discrete invariants of the topological insulator. 
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1. Introduction 



In condensed matter systems, most states of matter are classified by the symmetries they 
break. For example, a crystal breaks the translational symmetry, a magnet breaks the 
rotational symmetry and a superconductor breaks the gauge symmetry. These broken 
symmetry states are identified by the order parameter, and described by the effective 
field theory. The effective field theory is constrained by the broken symmetry, and it 
defines universality classes of different states of matter, and predicts universal physical 
properties. 

The quantum Hall (QH) state is a topological state of matter which does not fall 
into the conventional Landau paradigm. The long distance and low energy properties 
of the QH state is generally described by the effective topological field theory 
Ci 



S cS = — I d 2 xdte^ VT A ll d u A T , 

4:71 J 



(1) 



This topological field theory is generally valid for interacting systems, and describes the 
electromagnetic response in the long wave length limit [lj]. In the non-interacting limit, 
the Hall conductance a xy = G\ can be expressed as a topological invariant over the 
two-dimensional (2D) Brillouin zone[2j], given by 

Ci = ~ [dk x fdk y f xy {k) (2) 
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are the abelian Berry curvature and potential for the band states |ak). Time reversal 
(TR) symmetry breaking is essential for the QH effect. 

More recently, new topological insulator states have been theoretically predicted 
and experimentally observed in HgTe quantum wells, BiSb alloys, Bi2Te3 and Bi2Se3 
crystals^; fl; H; B 3 H ; 0]. These topological insulator states physically arises 
from the spin-orbit coupling in electronic structures, and are protected by the TR 
symmetry^; 0; [ll|; 0; [l3j]. A general theory describes the three dimensional (3D) 
topological insulator in terms of a effective topological field theory 12| , given by 
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where the last term is the topological term and a = e 2 /he is the fine structure constant. 
For a period system, the partition function is invariant under the shift of 9 by integer 
multiples of 2tt. Since TR transformation replaces 9 by —9, only two special, quantized 
values of 9 = and 9 = it are consistent with the TR symmetry. Therefore, all TR 
invariant insulators fall into the two distinct classes, disconnected from each other. The 



Equivalent topological invariants of topological insulators 



3 



Standard Model of elementary particles also admits such a topological term, therefore, 
the study of the topological insulator can also shed light on the fundamental topological 
interactions in Nature 14 . 

We can always couple an arbitrarily interacting electron system to the external 
electromagnetic field, and integrate out all the electronic degrees of freedom to obtain 
the effective field theory ([3]). For a TR invariant insulator, only two discrete possibilities 
of 9 = and 9 = tc can be realized. For a non-interacting band insulator, Qi, Hughes 
and Zhang (QHZ) gave an explicit formula for the 9 parameter as an integral over the 
3D Brillouin zone 

9 = 2nP 3 (9) = J d 3 ke^Tr{[^(k) - ^(k) • % (k)] • a k (k)} (4) 

where 

f Q/3 = da al3 - da a(3 + i\a- a -] af3 
af(k)= -*<a,k| A|0 jk ) 

are the non-abelian Berry curvature and potential for the band state k). This 
topological invariant has the physical interpretation of a magneto-electric polarization, 
which can be directly measured experimentally 12j; ll5 . 

We notice a beautiful symmetry between the topological field theory of the TR 
breaking QH state and the TR invariant topological insulator state. Both the space- 
time integrals (PQ), (j3j) and the Brillouin zone integrals (j2J), (jlj) are integral topological 
invariants in the theory of differential geometry. We see that the TR breaking QH state 
is described by a Chern-Simons integral over the 2 + 1 dimensional space-time (ED), and 
a first Chern integral over the 2D Brillouin zone (j2J) , while the TR invariant topological 
insulator is described by the second Chern integral over the 3 + 1 dimensional space-time 
([3]) and a Chern-Simons integral (HI) over the 3D Brillouin zone. These are the deepest 
and most natural topological invariants in mathematics, and it is gratifying to see that 
they also describe topological states realized in Nature. 

TR invariant insulators form an universality class extending over the 4D, 3D and 
2D space. The root state of this universality class is the topological insulator in 
4D 16|; Li?]] - In fact, it was the first TR invariant insulator state introduced theoretically, 
and historically it was referred to as the 4D QH state. It is described by a topological 
field theory 

ScS = ^I ^ xdteilupai "AAAAA T (5) 

Under TR transformation, 

A -> A , A, -> -At (6) 

therefore, we see that the topological field theory in 2 + 1 dimensions (OQ) breaks TR 
symmetry, whereas the topological field theory in 4 + 1 dimensions (JSJ) preserves the 
TR symmetry, and naturally describes the TR invariant topological insulators. For the 
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case of non-interacting fermions, the coefficient is given explicitly by the second Chern 
number 

f 
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The TR invariant topological insulator defined by (jSJ) and (jTJ) in the 4D space 
naturally generalizes the corresponding quantities of the TR breaking topological 
insulator defined by ([1]) and ([2]) in the 2D space. This is the reason why it was historically 
referred to as the 4D QH state [la; It is the root state for all TR invariant topological 
insulators in 3D and 2D, which can be obtained from the root state in 4D through the 
process of dimensional reduction 
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Top olog ical invariants in 3D band insulators have been studied from a different 
approach[10|; ll3j. In a beautiful series of papers, Fu, Kane and Mele (FKM) introduced 
a Z 2 topological invariant 0; Eol ; [l~8| for the strong topological insulator, expressed as 
a discrete product over the eight time- reversal invariant momenta (TRIM), explicitly 
written as 
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where 8{ is quantity defined at TRIM Tj, the detailed form of which will be given in the 
following sections. Since weak topological insulators may not be generally robust, we 
shall not be concerned with their definitions here. 

Therefore, there are now two topological invariants defined for the topological 
insulator. They are motivated by different logical reasoning and have different 
mathematical forms. The integral invariant (j4]) given by QHZ is physically measurable in 
terms of the magneto-electric polarization, and leads directly to the general topological 
field theory ([3]). The discrete invariant (jSJ) given by FKM has the distinct advantage 
that it can be easily evaluated, especially for cr yst als with inversion symmetry. Applied 



to concrete models of topological insulators 12j; |l9j], these two definitions yield the same 



result. However, it is highly desirable to prove the general equivalence between these 
two definitions. In this paper, we show explicitly that one can transform the integral 
invariant (j4]) to the discrete invariant (jSJ) exactly, proving the precise equivalence between 
these two definitions. 



2. A pedagogical example 

Because our approach involved the mathematical concept of degree of map, we shall 
give a brief introduction to it. Our presentation in this section is mainly designed for 
intuitive understanding, rather than mathematical rigor. Therefore, we shall illustrate 
the idea of degree of map in a simple one-dimensional example. For more mathematical 
details, c. /. Ref. [20J. Those readers who are already familiar with this subject may 
skip this section. 

Let us consider one concrete example. Consider a map / : M — > N, where M 
and N are both one-dimensional circle S 1 (see FigJT]). The coordinates of M and N are 



Equivalent topological invariants of topological insulators 



5 




denoted as <p and 9, respectively. The degree of map / is just the number of times that 
M covers N under / and it is often called "winding number" . 

The standard integral form of the winding number is defined as 

deg(/) = — / d0(<f>) = — -d<j ) = neZ (9) 

2lT J<f>=0 ZTT J<j>=o dq> 

which has a simple geometrical interpretation. For simplicity, let us assume that the 
point = maps to 9 = 0. When goes from to 2tt, 9 goes from to 2irn. This 
integral form of the winding number can be expressed in a discrete form. We arbitrarily 
choose a image point p on the image manifold N (see Fig{T]for illustration), and count 
the number of source points pi, p%, ps, etc on the source manifold M, where the source 
points are weighted by the +1 or —1 sign depending on the orientation of the map. For 
example, image point q has only one source point q, which maps onto q in a clockwise 
sense. Therefore deg(/) = 1. On the other hand, the image point p has three source 
points pi, p~2 and p%. pi and pi map onto p in a clockwise sense, while p~2 maps onto p in 
a counterclockwise sense, giving deg(/) = +1 — 1 + 1 = 1. We see from this example that 
the integral invariant of the winding number can be reduced to counting the weighted 
number of source points of a given image point - this is a discrete invariant for the 
winding number. 

The above picture can be generalized to higher dimensions. Generally, for a map 
/ : M —>■ N, where M and N are both d dimensional orientable manifold, we can 
defined the degree of map as 

deg(/) = / /» (10) 

where u is a d— form volume element on TV satisfying J n lv = 1, and f*{uj) is the pullback 
of io to M under /. In our one-dimensional example discussed earlier, u = d9/2n, and 
f*(uj) = d0(0)/27r = (d9/d<t>)d<t>/2ir. 

In our one dimensional example, we have seen that there is a discrete form of the 
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invariant deg(/). This can also be generalized to higher dimension. It is given as 20] 

deg(/) = N[f-\p), J f - Hp) > 0] - N[f-\p), J rHp) < 0] (11) 

where J is the Jacobian of the map, and N[f^ 1 (p), Jf-i(p) > (<)0] denotes the number 
of source points mapping to p with a positive (negative) Jacobian, respectively. Without 
going into the general derivation of this formula, we just give the explanation in our 
one dimensional example. In that example, Jf-^M is the direction sense (clockwise or 
counterclockwise) of the map at the point p, which is exactly what we have discussed 
earlier. We also mention that the point p in Eq-t fTTJ) should be regular, which means 
that the Jacobian of map are nonzero at Because the set of non-regular points 

has zero measure, we can always perturb the map to remove non-regularity at a given 
point. 

In summary, we showed that the degree of a map can be expressed in two equivalent 
forms, the integral form ( FlOi) and the discrete form ( TIT]) . In the next section, we will 
apply this idea to the integral and the discrete invariants of @ and ([8]) , which are degree 
of map deg(/) modulo 2, simply denoted as deg 2 (/). 

3. Equivalence between integral and discrete topological invariants 

Let us start from the band structure of time-reversal invariant topological insulators. 
To simplify the problem, we first assume that there are no degeneracies except those 
required by time reversal symmetry. Suppose that there are 2N filled bands labeled by 
(a, /?,••■), where the number of filled bands is even because of the Kramers degeneracy. 
The 2N by 2N matrix -B(k) is defined by 

| -k,a> =53S^(k)|e,k,/3> (12) 

P 

where |0,k, (3) = 6|k, (3) and 9 is the TR operator. The property 2 = —1 is crucial 
for TR invariant Fermi systems. By direct calculation we get the inner product 

(-k,e*|e,k,/3) = 5> a7 (k)<e,k, 7 |e,k,/?) 

7 
7 

= B aP (k) (13) 

and an important property of -B(k) 

B Ql3 (-k) = (k,a|e,-k,/3) 

= -<-k,/3|9,k,a) 

= -£/? a (k) (14) 

where the facts that 6 2 = —1 and (a\Q,j3) = — (/3|0, a) have been used. Therefore, 
B a p(— k) is anti-symmetric at the eight TR invariant momenta(TRIM), which enables 
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the definition of Pfaffian at these points. The definition of the discrete invariant given 
by FKM 
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181 ] is expressed as 
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with the Pfaffian of B written as Pf[Z?]. Although it seems that all the quantities 
appearing in this definition are local, global information on T 3 is encoded because a 
global basis of wavefunction is required in this definition. The existence of global basis 
is not obvious and we shall present a discussion on this in the appendix. The unitary 
matrix B(k) defines a map 

/ : T 3 -> U{2N) (15) 

This map can be simplified since we have the assumption of absence of accidental 
degeneracies. We divide the 2N filled bands into N TR pairs. Because energy eigenstates 
with different eigenenergy are orthogonal, we have B a p(k) — if | — k, a) and |k, f3) 
belong to different pairs. Therefore, all the inter-pairs elements of matrix B(k) are zero. 
Written explicitly, the matrix B(k) takes the following form 

B 1 (k) 

B 2 (k) 

£ 3 (k) 



B(k) 



(16) 



B N (k) 

where each B m (k) is an U(2) matrix. Therefore, the map / splits into N maps 

f m :T 3 ^U(2) m=l,2,---,N (17) 

The next step is reducing U(2) to SU(2), which is a three-dimensional manifold. 
The motivation of this step will be clear later when we consider the degree of certain map. 
To this end, we factorize B m (k) = exp[i6 l m (k)]-u m (k), where u m (k) is a SU(2) matrix 
function of k. Because of the facts that exp[i9 m (k)]u m (k) = — exp[i9 m (k)][—u m (k)} 
and — u m (k) is also SU(2) matrix, there is ambiguity in this factorization. To 
avoid this ambiguity, we can choose a point k in T 3 , and choose one factorization 
B rn (ko) = exp[i6 , m (k )]'U m (ko). The U(l) factors of other k points are determined by 



m (k) = m (ko) ~\ \ dkV k det[5 m (k)] (18) 

It is easy to check that with this equation u m (k) will have unitary determinant. To make 
the integral in the above equation unambiguous, we require that § t dkVkdet[5 m (k)] = 
for arbitrary loop I in T 3 . For contractable loops, this is trivially satisfied. For non- 
contractable loop, consider the non-contractable loop AOB in Figj2] as an example, we 
have 
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Figure 2. The loop (closed line) AOB in T 3 . The point O is k = (0, 0, 0). Points A and 
_B are the same point because the periodic boundary condition. AOB is homotopically 
nontrivial, i.e. AOB cannot be continuously deformed to a single point. 



I dkV k detLB m (k)] = / dkV k det[.B m (k)] + / dkV k detLB m (k)] 
J aob J ao Job 

= (19) 

which is a result of the relation det[B m (— k)] = det[S m (k)]. Therefore, we have a 
factorization of U(2) into SU(2) and U(l). With this factorization, we also have that 
$dkVk#m(k) = (for arbitrary loop I). Therefore, we can adiabatically deform the 
£7(1) factor to 1. Having got rid of the U(l) factor, the maps f m : T 3 — > U(2) are 
deformed to g m : T 3 — >• SU(2). This completes our discussion on matrix B(k). With 
the aid of -B(k), we can give geometrical interpretation of topological invariants of TR 
invariant topological insulators. 

In order to prove the equivalence, we first show how to relate the topological 
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invariant P 3 in Eq. (Tjy) to a winding number. Similar formula appeared in Ref. 
in a slightly different form but we would like to include the discussion here to make our 
discussion self-contained. We start from the integral topological invariants 

l r rr . ... 2 



167T 



J d 3 ke^Tr{[^-(k) - -iai(k) ■ % (k)] ■ a fc (k)} (20) 



Equivalent topological invariants of topological insulators 9 

which is valid independent of the TR symmetry. In a TR invariant system, with the aid 
of B(k), we have 

af(-k)= -^-k^lcUj-k,/^ 

= i x: B aa ,(Q, k, o'l^^ie, k, /?'» 

= i ^ B aa ,B* p/3l (0, k, a'l^J©, k, /?') + i ^ B aa ,d k .By6 a >p> 

a' /3' a'0' 

= ^B^K'^MYB^ + i^B^d^Blp (21) 

which can be simply written as 

Oi(-k) = B(k)a*(k)B^k) + ^(k^P^k) (22) 

Therefore, the field strength satisfies the following relation 

fij(-k) = d- ki aj(-k) - d- kj ai(-k) +i[ai(-k),aj(-k)] 

= - P(k)/*.(k)pt (k) (23) 

These results imply that the non-abelian Berry potential and curvature at — k are simply 
related to those at k by a non-abelian gauge transformation P(k). Therefore, we obtain 
that 

P3= vhJ rf3keU ' fcTr W- k ) - ^(-kK(-k)H(-k)} 

d 3 ke^Tr{[^(k) - ^(k) % (k)]a fc (k)}* 



16n 2 J -i WW v-/ 3 

^ | ^ke^Trc^Pa*^) 



" 2^2 / rf 3 ke^ fc Tr[( J B^ J Bt)( jB ^. jB t)( jB ^, jB t ) ] 
= -P s * - 24^/ rf 3 ke^Tr[(^ J Bt)(S9 J St)(S9 feJ Bt)] 

= - P 3 - ^ / rf 3 ke^ fc Tr[( J B,9, J Bt)( J B^ J Bt)( J B9 fcJ Bt)] ( 24 ) 

or 

2P 3 = / rf 3 ke^Tr[( J B9 iJ Bt)( jB ^. jB t)( jB4jB t)] (25 ) 

Both the LHS and RHS of the above equation depend on gauge choices, but the parity 
of 2P 3 is gauge invariant 

2P 3 (mod2) = -7^/ d 3 ke^Tr[(Pd i P t )(Pd ? P t )(Pd fe P t )] (mod2) (26) 

which is the important quantity for characterization of the topological insulator. 

In our previous discussion, we have split B(k) into direct sum of SU(2) 
matrices(P m , m — 1, 2, • • • , N). It is readily obtained that 

N 

2P 3 (mod2) = (mod 2) (27) 

m=l 
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where 

-2^S / rf 3 ke^TV[(S m ^St i )(S m a i St l )(S m a fe St i )] (28) 
We notice that v m is just the winding number, or the integral form of degree of the map 
g m : T 3 -> Si7(2) (29) 
The mod 2 degree of this map is given by 

deg 2 (5w) = ^ m (mod2) (30) 
Therefore we have 

N 

2P 3 (mod2)= £deg 2 (<? w ) (31) 

m=l 

This is the integral form of degree of map. The natural question is whether we can 
find a discrete form, as discussed in the previous section. Let us denote the image of 
k under the map g m as B m (k), which is a SU(2) matrix. As has been noted, B m (Ti) 
is anti-symmetric at TRIM Tj. There are only two anti-symmetric matrices in SU(2), 
given by 

M = I , ' ! • A, 





It is easy to see that Pf[^i] = 1 and Pf[v4 2 ] = —1. In principle, we can pick any image 
point on the 577(2) to perform the counting for the degree of the map. However, because 
of the TR symmetry, the counting is particularly simple if we pick either A 2 or A\. We 
now count the number of source points k, modulo 2, which map onto A 2 . We first notice 
that if B m (k) = A 2 , then B m (— k) = — -B^(k) = A 2 . Therefore, if k maps to A 2 , so 
does — k. If k is not one of the TRIM, k and — k are distinct, and these two points do 
not contribute to deg 2 (g m ), which counts the number of source points modulo 2. The 
only source points which may contribute to deg 2 (g m ) are TRIM, where k and — k are 
equivalent. It is also important to note that the calculation of mod 2 degree of map is 
easier than that of the integer degree, which involves the ± signs of the Jacobian of the 
map. Because — 1 = l(mod2), we do not need to be concerned about the signs of the 
Jacobian when calculating degree of map mod 2. 

Because we have reduced the map from U(2) to SU(2), det[.B(I\)] = 1, we have 
6 t = v /det[/i(r i )]/Pf[5(r i )] = Pf[B(Ti)] = U m ^[B m (Ti)}. Suppose that there are n m 
TRIMs which map onto A 2 under the map gj and 8 — n m TRIMs which map onto Ai. By 
counting the number of points which map to A 2 mod 2, we have deg 2 (g m ) = n m (mod2), 
and therefore 

(_i) 2p 3 = (-i)ELi de g 2 (^) = rj(-i)"™ (32) 



On the other hand, the discrete invariant defined in Ref. [10J is given by 

(-ir = unm)) = n pf [^( r *)] = ii(-^ nm m 
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Therefore, we proved the central result of this paper, namely the exact equivalence 
between the integral invariant of QHZ and the discrete invariant of FKM for the 3D 
topological insulator: 



In the proof given above, we made the assumption that no accidental degeneracy 
except the Kramers degeneracy occurs, which simplified the discussion. The extension to 
the generic case is straightforward because ir 3 (U(2N)) = n 3 (U(2)) holds for any integer 
N > 1. One can always deform the map / to one of the U{2) subgroups of U(2N), so 
that the proof discussed above applies. 
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Appendix A. Discussion on the global basis of wavefunctions 

In this appendix we shall show that the global basis of wavefunction exists on T 3 at 
presence of time-reversal symmetry. Because the global definition of B(k) and the 
discrete invariant Vq depends on the existence of such global basis, this point is important 
to our argument. 

We divide the 2N filled bands into N TR pairs. Within each pair, the Hilbert space 
is two dimensional at each momentum k. Therefore, the Hilbert space is naturally an 
U(2) fibre bundle on T 3 . Generically it is not evident that this bundle is trivial. The 
key to the existence of global basis is the TR symmetry. 

The "surface" of the cubic in FigfS] consists of three T 2 (two-dimensional tori). 
We consider the restriction of the U(2) bundle to one of these tori. The bundle on 
T 2 is trivial because the first chern number C\ = 0, which is a consequence of TR 
symmetry. Therefore, a global basis exists on the surface of the cubic. Such a basis can 
be extrapolated to the interior of the cubic since the the interior region is topologically 
trivial. We consider the transition function between the surface and the interior region 
on their overlapping region, which is topologically equivalent to S 2 (two-dimensional 
sphere). Because the second homotopy group 7r 2 (£7(2)) is trivial, mapping from S 2 to 
U{2) (the structure group of the fibre bundle) are all trivial. Thus the bundle is trivial 
on T 3 . Therefore, global basis does exist. 

For comparison, we note that in the case of integer quantum Hall states, the global 
basis does not exist on the Brillouin zone T 2 , because the relevant fibre bundle on T 2 
is non-trivial. From this example we also see that TR symmetry is necessary for our 
argument. 




(-1) 



(34) 
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